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We present a conditional experiment involving a para- 
metric amplifier and an avalanche photodetector to generate 
highly nonclassical states of the radiation field. The nonclas- 
sicality is robust against amplifier gain, detector efficiency 
and dark counts. At the output all the generalized Wigner 
functions have negative values, and this is exploited in or- 
der to reveal the nonclassicality through quantum homodyne 
tomography. 

Nonclassical states of light are relevant in many fields, 
which ranges from fundamental tests of quantum me- 
chanics to applications in quantum communication and 
measurements [Q. In the last two decades several 
schemes to generate different kinds of nonclassical light 
have been suggested, and some of them have been imple- 
mented . However, with the exception of squeezing, the 
generation and the detection of nonclassical light is ex- 
perimentally challenging. The aim of this letter is to sug- 
gest a simple experimental scheme to verify the quantum 
nature of light in a state as nonclassical as a Fock num- 
ber state. The main feature of our experimental scheme 
is its robustness against the possible imperfections of the 
setup, such as finite amplifier gain, nonunit detector effi- 
ciency and the occurrence of dark counts. 

Let us consider an active x^^^ crystal operating as a 
nondegenerate parametric amplifier (NOPA, for details 
of the experimental setup see the end of the paper). The 
NOPA, pumped at frequency up — uja + u>b, couples two 
modes a and b (idler and signal modes) via the medium 
nonlinearity. In the rotating-wave approximation, the 
evolution operator of the NOPA under phase-matching 
conditions can be written as U\ = exp [A (a^6^ — a6)] 
where the "gain" A is proportional to the interaction- 
time, the nonlinear susceptibility, and the pump inten- 
sity. For vacuum input, we have spontaneous parametric 
down-conversion and the output state is given by the 
twin-beam 



IV') - - E ^'»- ® l'^)^ ^ = tanhA. (1) 



n=0 



We now consider the situation in which one of the two 
beams (say, mode b) is revealed by an ideal avalanche 
on/off photodetector, i.e. a detector which has no out- 
put when no photon is detected and a fixed output when 
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one or more photons are detected. The action of an 
on/off detector is described by the two- value POVM 



Hi = I - Ho 



(2) 



being the quantum efficiency. The outcome " 1" (i.e 
registering a " click" corresponding to one or more incom- 
ing photons) occur with probability 



Pi = (T/'li^niii/.) = 
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and correspondingly, the conditional output states for the 
mode a is given by 
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In the Fock basis we have 
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The density matrix in Eq. (|^) describes a psewrfo-thermal 
state, where the vacuum component has been removed 
by the conditional measurement. Such a state is highly 
nonclassical, as also discussed in Ref. In the limit 
of low gain A <IC 1 the conditional state approaches 
the number state |1)(1| with one photon. The Wigner 
function W{a) of the state @) exhibits negative values 
for any gain A and quantum efficiency rjp^. In particular, 
in the origin of the phase space we have 
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One can see that also the generalized Wigner function for 
s-ordering Ws(a) = —2/{tts) J ^^71^0(7) exp[2/s|Q;-7p] 
shows negative values for s G (—1,0). In particular one 
has 



WM = - 
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If we take as a measure of nonclassicality the lowest index 
s* for which Wg is a well-behaved probability (regular, 
positive definite) Q Eq. (j^) says that for qi we have 
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s* = —1, i. e. gi describes a state as nonclassical as a 
number state. 

The fact that all the generalized Wigner functions have 
negative values may be exploited in order to reveal the 
nonclassicality of gi through quantum homodyne tomog- 
raphy. In fact, one has 
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and therefore |^ 



(8) 



(9) 



where (x) is the probability distribution of a random 
phase homodyne detection (with quantum efhciency r/n) 
and [Ws] (x) is the tomographic kernel for the opera- 
tor Ws , which is given by M 



Rr,AWs]ix) = 
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where $(a, 6; z) is the confluent hypergeometric function. 
Rjj^[Ws]{x) is a bounded function for s < 1 — r/^^ which 
represents the maximum index of the Wigner function 
H^s(O) that can be reconstructed by homodyne tomog- 
raphy with efficiency 7]^. In Fig. y we show a typical 
reconstruction of M^s(O) for different values of s, A and 
Va- 




FIG. 1. Reconstruction of the Wigner function in the ori- 
gin of the phase-space by (Monte Carlo simulated) homodyne 
tomography (sample of 5 x 10* data at random phase). Top: 
reconstruction of 14^-0.5(0) versus gain A, for r^A = 60% and 
r^H = 70%. Bottom: reconstruction of VF_o.3(0) versus the 
avalanche photodetector quantum efficiency t^a, for A = 0.5 
homodyne quantum efficiency rjH = 80%. In both plots the 
solid line is the theoretical value, and the inset shows the 
detection probability Pi. 

Besides quantum efficiency, i.e. lost photons, the per- 
formance of a realistic photodetector are also degraded 
by the presence of dark-count, i.e. by "clicks" that do 
not correspond to any incoming photon. In order to take 
into account both these effects we use the simple scheme 
depicted in Figure ||. A real photodetector is modeled 
as an ideal photodetector (unit quantum efficiency, no 



dark-count) preceded by a beam splitter (of transmissiv- 
ity equal to the quantum efficiency rj) whose second port 
is in an auxiliary excited state p, which can be a thermal 
state, or a phase-averaged coherent state, depending on 
the kind of background noise (a thermal or Poissonian). 
If the second port of the beam splitter is the vacuum 
D = |0)(0| we have no dark-count and the POVM of the 
photodetector reduces to that of Eq. (0). 
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FIG. 2. Model for a realistic photodetector. 

For the second port of the BS excited in a generic mix- 
ture i> = Y,s'^ss\s){s\ the P0VA4 for the on/off pho- 
todetection is given by (IIi = I — Bo) 



n + s 
s 



\n){n\ . (11) 



The density matrices of a thermal state and a phase- 
averaged coherent state (with M mean photons) are given 
by 
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In order to reproduce a background noise with mean pho- 
ton number N we consider the state i) with average pho- 
ton number Af = N/ (1 — 77a)- In these case we have 
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where T and P denotes thermal and Poissonian respec- 
tively, and Ln{x) is the Laguerre polynomial of order n. 
The corresponding detection probabilities are given by 
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For small N the two models lead to the same detection 
probability at the first order in N ~ Pf = Pi + 0[A2] 
with 



Pi = 
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In the following we will use the Poissonian background, 
and omit the index p. The conditional output state, after 
the observation of a click, is now given by 
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and the generalized Wigner function Wg (0) in the origin 
of the phase space is 
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The Wigner function W{0) for s = (as well as for s < 0) 
is no longer negative for all values of parameters. On the 
other hand, there is a large range of values of the quan- 
tum efficiency 77a and the gain A giving negative Ws{0) 
with s accessible by homodyne tomography with realis- 
tic values of the homodyne efhciency %. In other words, 
the generation and the detection of nonclassicality are 
robust also against the occurrence of dark counts in the 
avalanche conditioning photodetector. In Fig. ^ we show 
a typical reconstruction of ^^^(O) for different values of 
parameters 
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FIG. 3. Reconstruction of the Wigner function in the ori- 
gin of the phase-space by (Monte Carlo simulated) homo- 
dyne tomography (sample of 5 x 10^ data at random phase) 
Top: reconstruction of VK_o.3(0) versus gain A, for rj^ = 70%, 
rjii = 80%, and = 0.05. Bottom: reconstruction of VP'-o. 4(0) 
versus quantum efficiency r^A for A = 0.8, i]a ~ 75%, and 
— 0.08. In both plots the solid line is the theoretical value, 
and the inset shows the detection probability Pi. Increasing 
A'^ would generally need a larger threshold for tja 

In a practical implementation the NOPA consists of a 
type-II phase-matched KTP crystal that is pumped by 
the second harmonic of a Q-switched and mode-locked 
NdiYAG laser. Previously, such a NOPA has been em- 
ployed with parametric gains > 10 (|^|^ > 0.9), to gen- 
erate highly quantum-correlated twin beams of light at 
1064 nm [Q. By appropriately choosing the input quan- 
tum state, a similar setup was then used to demonstrate 



the production of squeezed-vacuum state with a high de- 
gree (5.8 ± 0.2 dB) of squeezing |^. In the present con- 
text, the twin beams, which are easily separable because 
of their orthogonal polarizations resulting from type-II 
phase matching, can be separately detected; beam a with 
a homodyne detector to verify nonclassicality, and beam 
b with an avalanche photodetector. The main challenge 
in the present experiment is the achievement of high de- 
grees of overlap (mode-matching efhciency) between the 
down-converted and the LO modes. Such overlap is non- 
trivial in pulsed, traveling-wave experiments owing to the 
distortion of the down-converted modes that is caused 
by the spatio-temporally Gaussian profile of the pump 
beam. With suitable choice of LOs, however, rjh > 70% 
has been achieved an adequate value for the present 
experiment (cf. Figs. ^ and In the measurements on 
beam b, the main challenge will be the selection of the 
appropriate mode, which can be performed by exploiting 
quantum-frequency conversion, a process that has been 
previously demonstrated ]lO[ |. 

In conclusion, we presented a robust experiment to ver- 
ify the quantum nature of light. This goal is achieved by 
a conditional scheme, where one of the entangled twin- 
beam exiting a NOPA is revealed by an avalanche pho- 
todetector, leaving the other one in a pseudo-thermal 
state with no vacuum component. The nonclassicality, 
as well as its verification by homodyne tomography, are 
robust against amplifier gain, detector efficiency and dark 
counts. 
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